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1. Introduction 



The study of the dynamics of driven interfaces in random materials is relevant to a 
wide class of physical problems, from fluid invasion in porous media to the motion 
of domain walls in random magnets [^0. The interface dynamics in these systems is 
controlled by the competition between an external field which exerts a driving force per 
unit area F on the interface and the pinning by impurities which impedes the motion. 
The interface is often modeled as an elastic medium that is distorted by disorder, but 
cannot break. The dynamics of such elastic interfaces driven through quenched disorder 
in the absence of thermal fluctuations has been studied extensively both analytically and 
numerically [p^p!0|. At zero temperature there is a sharp transition from a sliding state 
above a critical driving force Ft to a pinned state below Ft- The transition has been 
described as a critical phenomenon in terms of scaling laws and critical exponents |]9|,p!0 



The critical exponents depend on the dimensionality of the interface and of the embedding 
space, as well as on the geometry of the quenched disorder. Thermal fluctuations are 
expected to round the transition. Closely related systems that exhibit the same type of 
nonlinear collective dynamics and have been studied very extensively are charge density 
waves (CDW's) in low dimensional conductors [6-11]. 

The dynamics of weakly pinned flux lines in type-II superconductors is another prob- 
lem in this general class [^. A transport current density J flowing through a type-II su- 
perconductor in the mixed state in the plane normal to the external field exerts a Lorentz 
force per unit volume F ~ J on the vortex lines, which then moves across the current 
causing electric fields proportional to the vortex velocity and hence resistance. Impurities 
and other defects in the material act to pin the vortices and impede their motion. Point- 
like material impurities, such as O2 vacancies, yield a quenched random pinning potential 
with short-range correlations. 

At low temperatures and fields, when the characteristic pinning energy barriers for a 
single vortex line exceed the energy associated with intervortex interactions, dissipation is 
controlled by single vortex dynamics. The problem can then be modeled as that of a one- 
dimensional elastic interface (string) driven through a random medium. It is important 
to distinguish between two types of short-range disorder that are relevant in different 
physical systems. These are usually referred to as random-field (or random-force) disorder 
and random-bond (or random-potential) disorder. The disorder is of the random-field type 



for domain walls in disordered ferromagnets [|T6[ or fluid interfaces in porous media |T^] . In 
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this case the pinning energy is the sum of the contributions from aU the impurities in the 
area or volume spanned by the invading fluid during its motion. In contrast, the disorder 
that controls the dynamics of domain boundaries in ferromagnets with random exchange 
interactions or flux lines in type-II superconductors is of the random bond or potential 
type. In this case the pinning energy is determined by the instantaneous position of the 
interface and it arises only from the impurities in its vicinity. The static behavior, i.e., at 
finite temperature and zero driving force, has been studied extensively and is known to 
depend strongly on the type of disorder [TB,rS|. In contrast, Narayan and Fisher recently 



argued that the critical behavior of driven interfaces at zero temperature is essentially 



independent of the type of disorder [|T0[ . Numerical work has not yet addressed this point 
conclusively. 

In this paper we present results of numerical simulations of the dynamics of an elastic 
string in quenched disorder of the random bond type. Our earlier studies of this model 
at zero temperature have focused on the critical properties at the depinning transition 
and the dynamics in the sliding state [rn,5[. We found that at T = there is a sharp 



transition at a threshold force Ft from a state where the string is pinned for F < Ft to a 
sliding state for F > Ft- The mean velocity vanishes as the threshold is approached from 
the shding state as ~ with f = \F - Ft\/Ft and /3 = 0.24 ± 0.1 |2g. The data can 
also be fit by a form v ~ l/ln\f\. The spatial range of velocity correlations in the sliding 
state is determined by a correlation length ^ that diverges at threshold as ^ ~ f~'^, with 
u = 1.05 ± 0.1. Due to the effect of random forces the driven interface profile is rough and 
can be characterized by a roughness exponent defined by 



< 



u{z, t) — w(0, t) 



>~ \z 



1^^, fl.l) 



where u{z, t) is the instantaneous position of a point on the string. Our numerical work 
at T = yielded ( = 0.97 ± 0.05. Numerical studies of related lattice models have yielded 
( = 1.25 ± 0.01 [^. The continuous elastic model of interface dynamics was studied 



recently by Narayan and Fisher [10] at T = via a systematic expansion in e = 5 — (i, where 
d is the dimensionality of the embedding space and 5 the upper critical dimensionality. 
For the case d = 2 of interest here they found u = ( = 1 exactly to all order in e and 
(3 = 1 — 2e/9 + 0{e^) = l/3 + (9(e^). These results are consistent with those obtained from 
our numerical work. 

We have now extended our work in two ways. First we have considered the string 
dynamics in the presence of both quenched disorder and thermal fiuctuations. Secondly, 
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we have investigated the properties of the zero temperature pinned phase below threshold. 
As expected, thermal fluctuations round the depinning transition and replace it by a 
smooth crossover. For temperatures well below a characteristic depinning temperature 
Tdp ~ (AiiTi^p)^/'^, where A is the strength of the quenched disorder (of range Rp) and K 
the string elastic constant, the temperature dependence of the mean velocity at the zero 
temperature threshold Ft scales with temperature as v{Ft-,T) ~ T^/'^, where (3 is the 
T = depinning exponent and r is a nonuniversal exponent, with t = | for continuous 
pinning potentials (see Figure 1). This scaling form was predicted some time ago [|I4[ by 
D. S. Fisher through a mean field theory for the Fukuyama-Lee-Rice model of CDW's. It 
was also recently verified numerically for both continuous and discrete CDW models in 
two and three dimensions [|22| . It is not a priori obvious that interfaces with short-range 
pinning potentials should behave like models of CDW's, where the pinning potential is 
periodic - and therefore correlated - along the direction of motion. In fact the critical 
exponents describing the T = depinning transition are different in these two classes 
of models. On the other hand, our results indicate that in both cases the dynamics at 
finite temperatures occurs via thermally induced hopping over small energy barriers and 
is insensitive to the details of the pinning potential, provided the time scale of the thermal 
hop is sufficiently short. The mechanism for destabilization of the soft modes of the system 
is the same for models with short-range pinning potentials and CDW's and the mean field 
argument proposed by Fisher applies in both cases. The details of the pinning potential 
only enter through the T = depinning exponent, /3. The string dynamics for T > T^p 
is qualitatively different and the scaling just described fails. This difference is clarified by 
the high force perturbation theory discussed in Section 3. Numerical studies of interface 
dynamics in 1 + 1 dimensions at finite temperature were carried out recently by Kaper 
et al. Our work on thermal effects complements this earlier studies by focusing on 

the scaling of velocity with temperature for T < Tdp, that was not discussed in [^, and 
on the qualitative difference in the string dynamics at low (T < T^p) and high (T > Tdp) 
temperatures. 

We have also investigated the properties of the zero temperature pinned phase below 
threshold. The response of the string when the driving force is increased towards thresh- 
old from below is dominated by localized forward jumps or "avalanches". The number 
distribution of avalanche sizes S\\ in the string direction, D{S\\]f), was conjectured to 



obey a scaling form near threshold |T^, characterized by a diverging correlation length 

DiSff) = -^D{S\^/^-l (1.2) 

where -0(511; f)dS\\df is the number of avalanches of diameter between S\\ and S\\ + dS\\ 
that occur when the reduced force is changed from / to / + df. Also, D{x) is a scaling 
function that decays rapidly for x >> 1. Using an e expansion in 5 — e dimensions, Narayan 
and Fisher predicted k = and u_ = u = 1, with u the correlation length exponent when 
the transition is approached from above threshold [ pi][ . 

We have studied the scaling of avalanches numerically by evaluating the distribution 
©(S'li) of avalanche sizes obtained in response to a small increase of the driving force 
for all forces below Ft, 

P(5||)= r dfD{S\\-J)^-^, (1.3) 

J-1 ^11 

and we find k' = 1.0 ± 0.2. The exponents k and k' are related by 

k'^k + —. (1.4) 

We find k' = 1.0 ± 0.2 and k = 0.05 ± 0.05, corresponding to u_ = 1.13 ± 0.30, consistent 
with the theoretical values. The value k = indicates that avalanches of large size are 
likely and the description of the interface or as an elastic string breaks down. This is also 
consistent with our result that the roughening exponent ( defined as A ~ 'S'||^'', where 
A is the area of an avalanche, is ^ = 1-0 ± 0.05. As discussed by Coppersmith in 
the context of CDW, large strain may develop in the d = 2 case of interest here. In a 
real system the string will relax these strains by pinching off a vortex loop around an 
impurity, a mechanism that is excluded in our model. Overhangs and loops could modify 
qualitatively both the dynamics in the pinned region and the T — depinning transition. 
Further studied are needed to address their role. 

The remainder of the paper is organized as follows. After briefiy describing the model 
in section 2, we discuss simple dimensional estimates and the results of the high field 
perturbation theory in section 3. The numerical result at finite temperature are presented 
in section 4, while the scaling of avalanches in the T = pinned state is discussed in section 
5. 
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2. The Model 

The specific model we have studied is an elastic string embedded in two dimensions 
and on the average aligned with the z direction (in superconductors this is the direction 
of the applied magnetic field). We describe the interface by its displacement u{z,t) in the 
direction of the applied force per unit length, F. The displacement is assumed to be a 
single- valued function of z at any time t, i.e., we ignore overhangs in the string and vortex 
loops that could be pinched off from it during the motion. In the absence of driving force 
the Hamiltonian of the string is the sum of the elastic energy of the string and a pinning 
potential from random impurities in the medium, 

where L is the size of the system in the z direction and K the elastic constant. Assuming 
purely relaxational dynamics, the equation of motion for the string in the presence of a 
constant driving force F per unit length in the transverse direction is given by 

5H 

ld^u{z,t) = --+F + rjiz,t) 

=Kdlu{z, t) + Fp{u; z) + F + r]{z, t), 

where 7 is a friction coefficient, Fp{u;z) = —dV{u; z)/du is the pinning force per unit 
length, and r]{z, t) is a Gaussian correlated Langevin force per unit length describing 
thermal noise, with < r]{z, t) >= and correlations 

< r]{z, t)r]{z', t') >= 2-iT5{z - z')5{t - t'). (2.3) 

The angular brackets < • • • > denote the thermal average. In the numerical model the 
pinning potential is written explicitly in terms of the interaction of the string with Np 

— * 

short-ranged pinning centers randomly distributed at positions Ri in the plane, 

V{u-z)^Y.^{\t-R,\), (2.4) 
1=1 

where r = {u{z, t),z). The interaction U (|r — Ri\) of the string with the i*'^ pin is approxi- 

— * 

mated by a potential well centered at the pin location Ri of finite range Rp and maximum 
depth Uq. The pins are uniformly distributed with areal density Up. When the mean pin 
spacing is large compared to the potential range Rp, the pinning potential can also 
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be described as a continuum gaussian random variable, with mean and correlations given 
by, 

V{u- z)V{u'- z') ~ AS{u - u')S{z - z'), (2.5) 

where 

A ^ {UoRp)^npRl{l + OiupRl)} (2.6) 

and the overbar denotes an average over disorder realizations. 

The overall motion of the string is described by a "center of mass" velocity, defined 

as ^ 

Vcn^{F,t) = y [ dzv{z,t), (2.7) 
^ Jo 

where v{z, t) = dtu{z, t) is the instantaneous velocity of a point on the string. The center 
of mass velocity is a fluctuating quantity since it depends of both the random positions of 
pins and the thermal noise. The average or drift velocity of the string is given by 

v{F) = < VcUF,t) >. (2.8) 

In the numerical calculation the average over different realizations of disorder (denoted by 
the overbar) is performed by averaging over time, since as time evolves the string samples 
different impurity configurations. 



3. Dimensional estimates and perturbation theory 

Considerable insight on the behavior of the driven string can be gained by simple 
dimensional estimates and by a perturbation theory in the pinning force. Most of the 
results discussed in this section have been obtained elsewhere for a general d-dimensional 
interface |2^. It is, however, instructive to summarize them here for the 1 + 1 dimensional 
case of interest. 

The relative importance of thermal fluctuations and quenched disorder in governing 
the equilibrium properties of elastic interfaces, in the absence of external drive, has been 
studied extensively. Thermal fluctuations always dominate at small length scales. In 
this regime they are responsible for small amplitude vibrations of the interface within a 
given pinning well and therefore lead to a smoothing of the pinning potential. At large 
length scales the importance of thermal fluctuations depends on dimensionality. For a one 
dimensional interface in d = 2 the disorder is always dominant at large length scales. This 
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is because the value of the roughness exponent C,eq of an elastic string disordered only by 
random bond quenched disorder, C,eq{d = 2) =2/3, is larger than the corresponding value 
for a thermal string in equilibrium, = 1/2. 

A disordered elastic interface in the absence of thermal fluctuations is characterized 
by the Larkin-Ovchinnikov (LO) collective pinning length, Lc- This is estimated by dimen- 
sional analysis by considering the energy fluctuation associated with displacing a segment 
of string of length L by a transverse distance Rp across a single pinning well, 

5F = + {UoRp)^^ih:R^. (3.1) 



The LO pinning length is obtained by minimizing (|3.1| ), with the result 

1/3 



Rpi^K/Uo-^UpRp) . When Lc » Rp, the string is pinned collectively by many pins. 
It is this weak pinning regime that is relevant for flux lines in superconductors. When 
thermal fluctuations are important the mean thermal displacement of the string exceeds 
Rp and the string experiences a random potential averaged over the length of its root mean 
square thermal excursion about equilibrium, < u^^ >^/^. The temperature dependent pin- 
ning length is estimated by considering the energy associated with a fluctuation of length 



L and transverse distance < >^/^, 



,2 



SF = + UoRp 



1/2 



(3.2) 



where < > is determined by requiring K < > /L ^ T. We flnd Lc{T) ^ 
LciT /Tdp)^ , where the depinning temperature T^p is deflned by < u^^ >^ Lc{T)T/K = Rp 
, which yields T^p ~ {UoRp)'^/^{npRpK)^/^. The temperature-dependent pinning length is 
deflned by interpolating between the T = LO length and the flnite temperature result, 

L,{T) = L41 + {T/Tdp)% (3.3) 

A more rigorous derivation of these results can be found in ||2^ . 

Considerable insight on the competing roles of quenched disorder and thermal noise 
can be gained by considering the string dynamics at large velocities, where the disorder can 
be treated as weak. For F » Ft , the effect of pinning is negligible and the string advances 
uniformly, with v ^ vq = F/'y. The deviations from this asymptotic behavior can be 
studied by a perturbation theory in F/ Fp that was introduced flrst by Schmid and Hauger 



and by Larkin and Ovchinnikov for T = and recently extended by Vinokur et al. [p6| , p7 | 
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to incorporate thermal fluctuations. To carry out the perturbation theory it is convenient to 
consider a frame of reference moving at the mean velocity v of the string. The instantaneous 
position of a point on the string is then written as u{z, t) = vt + Up{z, t) + uth{z, t), where 
Up{z,t) is the deformation due to quenched disorder, treated as small, and uth{z,t) is the 
contribution from thermal fluctuations to the displacement from the uniform sliding state. 
It is defined as the solution of, 

dtuth{z,t) = Kdluth{z,t) + r]{z,t). (3.4) 

The mean velocity is then given by 

1 



5v = - < Fp{vt + Up + uth\ z) >, (3.5) 
7 

with 5v = V — vq and 

-fdtUp{z,t) = Kdlup{z,t) + Fp[vt + Up{z,t) + uth{z,t)). (3.6) 



The right hand side of Eq. (|3.5| ) is evaluated in perturbation theory by treating Up as small. 
The first order correction to the average string velocity as compared to the asymptotic value 
f = f is given by the self-consistent solution of the equation 



— = -/ / G{0,t)Ss{k,t), (3.7) 



where p{k) is the single-vortex form factor that provides a short-length scale cutoff at 
k ~ ^/Rp, with k the wavevector in the direction of motion. Also, G{0,t) is the vortex 
Green's function evaluated at ^ = 0, given by 

G{0,t) = e{t) [ ^le-5'-^V7. (3.8) 
J 27r7 

The prime on the integration over the wavevector — oo < q < +oo in the z direction 
denotes a short distance cutoff at 1^1 ~ 2n/Rp which controls the crossover to single 
particle behavior at short times. Thermal fluctuations are responsible for the appearance 
in Eq. (|3.7| ) dynamical structure factor Ss{k, t) of a single vortex in the absence of disorder, 
given by 



In the absence of thermal fluctuations Ss{k,t) = 1 and the cutoff kg on the /c-integration 
in Eq. ( |377| ) is provided by the form factor p{k), i.e. ko ~ l/Rp. At finite temperature this 
cutoff is replaced by ko ~ [Rl+ < u^hi^) > /2]~^/^. 

By using the fiuctuation-dissipation theorem the mean square thermal displacement 
can be expressed in terms of the Green's function given in Eq. (|3.8D , with the result 

,2 m [' ^l^^-.-K^th 

K J 27rQ2 



<<,it)>=— -L [l-e-^'^y^ . (3.10) 



The (j-integration on the right hand side of Eq. ( |3.10|) is easily carried out, with the result 

T 

< ulhit) > ^ -y^Vt, t » to, 

4^ ^ (3.11) 
^ ——t, t « to, 

where to = R^/Dq, with Do = K/'j the diffusion constant of a point vortex, is the time 
scale for diffusion across the range Rp of the pinning potential. For t « to different bits 
of string of length < Rp are essentially uncorrelated along the z direction and one recovers 
single-particle behavior, < u^f^it) >~ t. 

As discussed earlier, at large times or long length scales disorder is always important 
and the long time divergence of Eq. ( p.ll|) is cut off by disorder. The effect of disorder 
can be approximately incorporated in Eq. ( 3.10 ) by introducing a large distance cutoff at 



q ~ 2n/Lc{T) in the q-integration. Neglecting for simplicity the crossover to single particle 
behavior, we obtain three distinct regimes 
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< w^Jt) > <Rl, to<t<tph = to {RpK/T 

< u1^{t) > ^ ^y^, tpH<t<td = tpH(T/Tdpy\ (3.12) 

We have introduced here two new time scales. The time tph characterizes the time scale for 
small phonon-like vibrations of the string within the potential well of a single pin, while ta 
is the time scale where disorder becomes dominant. The intermediate regime in Eq. ( |3.12| ) 
only occurs provided tph < td which corresponds to T > Tdp. 

We now proceed to approximately evaluate the correction 5v given in Eq. (|3.7| ) . The 
right hand side of Eq. ( |3.7| ) can be reduced to a one-dimensional integral that can be 
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/+00 ^j, fl/kv 
— k^ / dtv^e-'='[«p+<"-W>/2l. (3.13) 
-oo 27r Jq 



evaluated numerically. It is, however, more instructive to simply carry out a dimensional 
estimate of the integral. Neglecting for now the short distance cutoff in the g-equation in 
Eq. ( |378| ) , we obtain ^(0, t) = Q{t)^y^yJnKt. Inserting this on the right hand side of Eq. 
( p.7| ) and noting that the main contribution to the time integration comes from kvt ~ 1, 
we obtain 

5v _ A /•+°° dk , , /•i/'^" 

As the driving force F - and therefore the string mean velocity v - increases, the time cutoff 
Xjkv decreases. Inserting on the right hand side of eq. ( |3.13|) the form for < u'^f^{t) > from 
Eq. ( p.l2|) appropriate to each time regime and carrying out the integration, we obtain 

Sv A 1 Rp 

Sv 2A f 1 , ^ Rp 



V 15 \ nT''' j v^/^ ' V 

7/2 



tph< — < td, (3.14) 



Sv 2A 1 f e \ 1 . R 



or 



T2 K 

<F < 



1 - 


V 


^ F' 


-3/2 




vo 






1 - 


V 


F' 


-1/3 




vo 






1 - 


V 


F' 


-3/2 











KRf, Rp' 



T3 /T^ \10 



F < 



KRl 



For F > K/Rp, corresponding to the case where the time cutoff is realized at times shorter 
than to, the string is sliding so fast that correlations among different string elements have 
no time to develop and one recovers single particle behavior, with 6v/F ~ F~^. The 
results summarized in Eqs. ( |3.15| ) show that the shape of the v-F curve is qualitatively 
different at high and low temperatures. The perturbation theory yields f /F ~ 1 — CF~"^, 
where C is a constant that depend on temperature and the value of a in the various 
regimes can be inferred from Eqs. (|3.15 ). The curvature of the v-F curve is determined 



by ^r§^ ~ —a{a — 1)F~"~^ and is positive if a < 1 and negative if a > 1. If T << Tdp 
the intermediate region described by the second of Eqs. (|3.15|) does not occur and the 
perturbation theory yields ct = 3/2 in the entire region where perturbation theory applies, 
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up to F ~ K/Rp. The v-F curve has a negative curvature in this regime and thermal effects 
only affect the coefficient of the correction 5v. When T > T^p there is an intermediate 
region described by the second of Eqs. (|3.15| ) where a = 1/3 and the v-F curve has 
a positive curvature. This behavior is apparent in our data described in the following 
section(see Fig. 2). 

4. Numerical Results at Finite Temperature 

We have integrated numerically the discretized version of the equation of motion 
( p.2| ) for a string composed of discrete elements, each of dimensionless size Rp in the z 
direction. Here Rp is chosen as the unit of lengths. All forces are measured in units of the 
string tension K. The string elements interact via nearest neighbor elastic forces and are 
constrained to move only in the direction of the driving force. This model is appropriate 
for flux lines in layered superconductors, where each flux line can be thought of as a stack 
of interacting two-dimensional "pancake" vortices residing in the Cu02 planes. Periodic 
boundary conditions are imposed in the z direction. Each string element is subject to 
attractive potential wells of finite size Rp = 1 and maximum depth Uq in both the z and 
u directions, centered at the randomly distributed pin locations. The displacements of the 
discrete string elements are treated as continuous variables and the coupled equations of 
motion are integrated by a fourth order Runge-Kutta algorithm with a time step much 
smaller than the typical time to cross a single potential well (typically At ~ 0.1, where 
time is measured in units of to = Rp/Do). As in our T = simulations, we have been 
able to obtain reliable data at very small velocities thanks to our "pinning cells" method. 
The method consists in dividing the plane in "pinning cells" of dimension Rp. At T = 
the string can only move forward and each section of the string only "sees" the disorder in 
the pinning cell that neighbors it in the forward direction of motion. Thermal fluctuation 
can also kick the string to move backward, in the direction opposite to that of the driving 
force. At every iteration each string element needs to know the disorder in the cells that 
it has left behind, since these may be revisited. We have developed an algorithm to store 
at each step an adjustable number of "past" pinning cells for each section of the string, 
in addition to the "future" cell. The list of past cells is updated in parallel at each time 
step. For driving forces within 10% of the zero temperature threshold, the length of 
the simulation usually exceeds 10^ time steps, while shorter simulation times gave good 
averages for forces further from threshold. We have investigated systems of size from 256 
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to 16384. The computations were performed on the Connection Machine CM-2 and CM-5. 
In our dimensionless units the parameters of the model are the dimensionless pinning force 
Fp = Uq/K, the dimensionless areal density of pins, p = UpRp, and the dimensionless 
driving force F = FRp/K. In these units the temperature is measured in units of KRp 
and the depinning temperature is given by Tdp = Tdp/K = py^p^/^ and the LO collective 
pinning length is Lc = Fp yielding a dimensionless estimate for threshold force 

Ft = pFp and Ft = p^^^Fp^^ for strong and weak pinning, respectively. We will always 
refer to the dimensionless quantities below and drop the tilde to simplify the notation. 

Our simulations were carried out for Fp = 1 and p = 0.1, yielding Lc = 1.5 and 
Tdp = 0.46. The dimensional estimate for the threshold force gives Ft = 0.1, which is 
consistent with our result Ft = 0.2435 ± 0.005 obtained from the simulations. The mean 
velocity of the string is shown in Fig. 2 for various temperatures and system sizes. For very 
large driving forces the effect of pinning is negligible and the string advances uniformly, 
with V ^ F. For T ^ there is no sharp transition and the velocity at low driving forces 
is small but finite. The velocity versus driving force (v-F) curve always exhibits a tail 
with positive curvature below the T = threshold. At higher driving forces the curvature 
of the v-F curve changes sign and eventually approaches the asymptotic limit v ~ F. At 
T = 1.0 X 10-2, well below our dimensional estimate of T^p ~ 0.46, thermal effects have 
already washed out completely the depinning transition. Thermal effects are quantitatively 
more pronounced than expected from dimensional estimates. 

We distinguish three regimes characterizing the string response at finite temperature. 
At high driving forces the deviation from the asymptotic behavior v ~ F are well described 
by the perturbation theory discussed in section 3. If T << Tdp our simulation agree with 
6v/F ^ F~'^/2(with a crossover to 6v/F ~ F~'^ at very large driving force) and the v-F 
curve has a negative curvature in this regime. At higher temperature the v-F curve shows 
an intermediate region with positive curvature where Sv/F ~ The qualitative 

difference in the shape of the v-F curve at high and low temperature is consistent with the 
results of the perturbation theory discussed at the end of Section 3. The value of T above 
which thermal effects change qualitatively the response is, however, much lower than the 
value obtained for Tdp from dimensional estimates. Finally, at very high forces, 1 — v/F ~ 
F~2. This crossover is controlled by the range of the pinning potential correlations in the 
z direction. 

For driving forces within a few percent of the zero temperature threshold and suf- 
ficiently low temperature, the mean velocity exhibits the scaling behavior proposed by 
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Fisher in the context of CDW's. He argued that thermaUy-induced "hops" over pinning 
energy barriers are analogue to jumps resulting from ramping up the force at T = 0. 
The velocity evaluated at the T = threshold is found to scale with temperature as 
v{F = Ft, T) ~ T^/'^, with /3 = 0.24 and r = 3/2, as shown in Fig. 1. In addition in the 
region of low T and F close to Ft, our data can be fit to the scaling form proposed by 
Fisher 

z;(/,T)~T^/-S(/T-V-), (4.1) 

as shown in Fig. 3. Here B{x) is a scaling function that behaves as B{x) ~ for x oo. 
This scaling fails, however, at higher temperatures. 

Finally, at very small driving forces (F << Ft) and low temperatures (T << Tdp) the 
dynamics occurs via creep over pinning energy barriers and v ~ e"^*-^-*/^. The pinning 
energy barrier U{F) has been predicted to diverge when F ^ as U{F) ~ F~^ . The creep 
exponent n has been estimated by dimensional analysis as = 1/4. It should, however, 
be noticed that this result was obtained by assuming that the roughening exponent of the 
string has the equilibrium value C,f,q = 2/3. Our numerical results in this region are not 
inconsistent with a small value of /U, as shown in Fig. 4, but are insufficient to either 
confirm or discard the assumption that the creep dynamics is controlled by the single 
diverging energy scale U{F) and to determine the creep exponent conclusively. A different 
approach may be needed to address this point. 

5. Avalanches in the T = pinned state 

In this section we discuss the critical behavior of the interface at zero temperature in 
the pinned region, as the threshold Ft is approached from below by slowly increasing the 
driving force. As discussed recently by Narayan and Fisher |jT^ for the interface problem 
and in more detail by Narayan and Middleton for CDW's, local instabilities occurs as 
F is increased towards threshold resulting in "avalanches" of various size. The notion of 
avalanches was introduced to describe the large response to a local perturbation in models 
exhibiting self-organized criticality (SOC)[^. Here we discuss avalanches obtained in 
response to a small increase of the driving force, which provides a global perturbation that 
affects equally all the discrete string elements. This perturbation, if sufficiently small, will, 
however, only trigger instabilities locally. The resulting response consists of discontinuous 
local jumps of portions of the interface (i.e., avalanches), with a distribution of sizes. The 
size of an avalanche can be defined in terms of its area (or moment)^ or of its diameter 
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i^ii in the direction of interface. We start the system in a pinned configuration at a driving 
force F < Ft- Increasing the reduced force / = F/Ft — 1 from f to f + df triggers a 
forward jump of sections of the string, until a new metastable pinned state is reached. 
The avalanche is characterized by the total area A swept by the string as a result of the 
increase in driving force, 

A = [ dzu{z). (5.1) 



Assuming u{z) ~ z^, with C a roughening exponent, the area and diameter of the avalanche 
are related by 

A ~ (5.2) 

Figure 5 shows a plot of the area of the avalanches versus their diameter. The straight line 
has slope 2 and we find C = 1 ± 0.05. 

Fisher and Narayan have conjectured that the distribution of avalanche sizes near 
threshold obeys the scaling form given in Eq. ( |1 . 2| ) ||1U|| . Using Eq. (|5.2|) , one can also im- 
mediately obtain the scaling form for the number distribution of avalanche areas, Da{A] /), 
given by, 

Da{A; f) = -j^DiA'/^'+^^/C-), (5.3) 



with 



= (5.4) 



Rather than considering the distribution of avalanches at a fixed force /, we have evaluated 
numerically the distribution of avalanche areas and diameters integrating over all driving 
force below Ft- The total distribution of avalanche size is given in Eq. ( |1.3| ). The 
corresponding distribution of avalanche areas is given by 

^a{A) = 1° dfDAiA; f) ~ (5.5) 

The corresponding exponents are related to the exponents defined in ([1.2|)and ( p.3| )by 

^ 1 

K — K ^ , 

1 (5-6) 



Ka = KA-\- 



^-(1 + 



To generate avalanches, we start with the string in a metastable pinned configuration 
well below Ft and study the response to a small increase AF of the force (we have used 
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AF = 1 X 10~^ in most of our simulations). The increase in force triggers local avalanches 
which are recorded, until the string reaches a new pinned configurations. The procedure 
is then repeated by stepping up again the force of an amount AF until Ft is reached. 
If, however, the string slides as a whole in response to the small perturbation, we discard 
this event, return to the original pinned configuration and repeat the procedure with a 
smaller AF. These "giant avalanches" are discarded because they are a finite-size effect, 
characteristic of the response of finite systems above threshold. The distributions are 
shown in Figs. 6a and 6b for two system sizes. They exhibit a power-law decay with 
k' = 1.0 ± 0.2 and k'^ = 0.5 ± 0.15. 

In order to gain further insight in the "shape" of the avalanches, we define the number 
distribution Du{A.u; /) of the local displacement advances Au following the ramping of the 
force. More precisely Du{Au; f)d{Au)df is the number of displacement advances between 
Au and Au -f d{Au) that occur when the reduced force is increased from / to / + df . To 
develop a scaling ansatz for Du{Au;z) we proceed as follows. Let n{Au\A)d{Au) be the 
conditional number of displacements between Au and Au -|- d{Au) that occurs within a 
given avalanche of fixed area A. We assume that the transverse shape of the avalanche (in 
the direction of motion) can be characterized by a single length scale defined as 

S^=A/S\\. (5.7) 

then conjecture a scaling ansatz for n{Au\A) of the form 

S 

n{Au\A) = -±h{Au/S\\), (5.8) 

where fi{x) is a scaling function that depends on the detailed shape of the string and the 
prefactor in Eq. ( p.8| ) is determined by the normalization condition 

d{Au)n{Au\A) = S\\. (5.9) 







The distribution of displacements -Du(Au; /) can then be written as 

Du{Au-f) = [ dADA{A;f)n{Au\A). (5.10) 



By inserting Eq. (|5.8|)in Eq. (|5.10| )and making use of Eq. (pT2|)and (p77|), we obtain 



D^{Au-J) (5.11) 
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with 

Hu = + ka{1 + ^). (5.12) 

Finally, the number distribution of the local displacements Au integrated over all forces 
up to Ft decays as 

^ A^' ^^-^^^ 

with 

< = «u + (5.14) 

The distribution Vu{Au) can be evaluated with excellent statistics and is shown in Fig. 7. 
We find < = 0.0 ±0.10. 

The primed exponents k', and k'^ governing the scaling of the distribution of 
avalanches for all forces below Ft are related by the same relations ( p.4| )and ( |5.12| )that 
hold among the unprimed exponents, 



''^^i (5.15) 

f^u = i^'a + ^{i^'a - !)• 

The exponents obtained from our numerics satisfy well these scaling relations. 

The correlation length exponent v_ can be inferred from ( p.6| ) or ( p.l4| ) if at least 
one of the distribution at fixed driving force is computed numerically. This is in general 
more difficult since the distribution is quite sensitive to the value of the force increment 
used, which needs to be made very small. We have evaluated numerically D{S]^\] f) for 
three different driving force. The distribution for / = 0.063 is shown in Fig. 8. We find 
K = 0.05 ± 0.05. Using the first of Eq. ( |5.6| ) and the value k' = 1.0 ± 0.2 quoted earlier, we 
obtain z/_ = 1.13±0.30. The result is again consistent with theoretical predictions. Finally, 
as a consistency check, we can now use v_ = 1.13 ±0.30 in the second of Eqs. ( p.6|) and Eq. 
( ^.14|) to obtain = 0.07 ± 0.22 and Ku = —0.95 ± 0.15. These in turn satisfy the scaling 
relationships (p.4|)and (|5.12| ) . 

The distribution of avalanches in discrete interface growth models has been studied 
numerically by Sneppen [|2^ and by Sneppen and Jensen [^^. In their model the growing 
interface is maintained in a "critical state" by a local growth rule similar to that used 
in invasion percolation models that prevents the formation of overhangs. The spatial 
and temporal correlations between successive growth events and avalanches in this model 
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have been studied extensively by Leschhorn and Tang |31]. They find a rather complex 
behavior where the distribution of growth events shows dynamical scaling only locally. This 
is because in the Sneppen model the driving force is self-tuned to maintain the interface 
in a "critical state" at the onset of steady state motion, thereby introducing additional 
spatial and temporal inhomogeneities in the model. 

In contrast, in our continuous model the distributions of various measures of avalanche 
sizes displays a well defined dynamical scaling characterized by a single correlation length 
^_ ~ |/|"~ as the threshold is approached from below. The correlation length exponent 
zv_ equals the correlation length exponent v obtained when the transition is approached 
from above. On the other hand, the values of the scaling exponent of the avalanche size 
distribution, k 0, and of the roughness exponent, C ~ I7 indicate that the interface may 
develop large local gradients, which were explicitly excluded in our model. A physical 
interface will relax these large strains by pinching off loops around the pinning centers, a 
mechanism that may modify qualitatively the dynamics in the critical region both below 
and above threshold. In particular this could reduced considerably the value of the rough- 
ness exponent Further studies of a more general model that incorporates strong elastic 
nonlinearities and allows for overhangs and loop generation are clearly needed to address 
these questions. 
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Figure Captions 



Fig. 1. Creep velocity as a function of dimensionless temperature T at the T = thresh- 
old force Ft- The sample size is indicated. The line shows the slope /3/t = 0.16 
for the exponents values discussed in the text. 

Fig. 2. Creep velocity as a function of dimensionless driving force at various tempera- 
tures. The straight line shows the asymptotic uniform velocity v ^ F. The data 
shown summarize results for three different system size : L — 1024, 4096, 16384. 
The same symbol is used for different system size as no significant size dependence 
was observed. 

Fig. 3. Scaled creep velocity vT~^/'^ versus scaled reduced force fT~^^'^ for the three 
lowest temperatures of Fig. 2, using r = 1.5 and (3 = 0.24. 

Fig. 4. Log-log plot of the energy barrier U{F) defined as U{F) — Tlog{F/v) versus 
driving force in the low temperature creep region(from T = 3.16 x 10~* to T = 
1.00 X 10-2). The dashed line has slope -1/4. 

Fig. 5. The area of avalanches obtained by successively stepping up the driving force by 
small increments / ~ 4.0 x 10~^ for — 1 < / < versus their diameter in the 
direction of the interface for both L = 1024 and L = 4096. The straight line has 
slope 2. 

Fig. 6. The distribution of avalanche diameters (a) and areas (b) integrated over all 
driving forces. The straight line has slope -2 for (a) and — | for (b). 

Fig. 7. The number distribution of displacement increments Au for all driving forces 
below threshold for L = 4096. The straight line has slope -1. 

Fig. 8. The number distribution of avalanche diameters for fixed driving force / = 
-0.063. The dimensionless AF used here is AF f=i 7.0 x 10"^. The straight 
line has slope —1.05. 
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